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SUMMARY 
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Analytic  equations  were  determined  for  the  force  and  moment  coeffi¬ 
cients  and  stability  derivatives  of  various  shapes  in  Newtonian  flow.  These 
results  follow  from  an  application  of  the  general  method  presented  in  the  "First 
Quarterly  Technical  Report". 

Equations  for  all  the  important  damping  derivatives  of  the  cone  at 
zero  sideslip  and  roll  have  been  determined  as  a  function  of  angle  of  attack. 

The  constants  in  these  simple  equations  can  be  determined  numerically  for  any 
C.G.  location  and  any  cone  angle.  For  slender  cones,  with  Certain  restrictions 
on  the  C.G.  location,  the  constants  in  the  stability  derivative  equations  have 
also  been  determined  analytically.  Similar  results  have  been  obtained  for  the 
delta  planform  with  diamond  shape  cross-section  perpendicular  to  the  longitudinal 
axis.  A  preliminary  investigation  has  also  been  conducted  on  the  elliptic  cone, 
and  equations  have  beer,  established  for  the  longitudinal  force  coefficients. 

A  preliminary  investigation  has  been  made  of  the  effects  of  bluntness 
on  the  constants  in  the  stability  equations.  These  results  have  been  applied 
to  the  determination  of  analytic  equations  for  the  longitudinal  force  coeffi¬ 
cients  of  the  blunt  cone. 

An  initial  investigation  has  been  undertaken  into  other  hypersonic 
theories,  besides  Newtonian,  for  predicting  surface  pressures.  This  investi¬ 
gation  is  continuing.  Based  on  these  preliminary  results,  a  method  has  been 
outlined  for  determining  the  variations  of  the  hypersonic  stability  derivatives 
with  Mach  number.  The  method  will  be  checked  against  experimental  results 
during  the  next  quarter. 
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All  symbols  riot  shown  here  are  defined  in  the  "First  Quarterly  Technical  Report" 

Ag  Base  area  of  the  bottom  half  of  a  spherical  cap  segment. 

C  Pressure  coefficient 

P 

e  Ratio  of  horizontal  to  vertical  axes  of  a  cone  of  elliptic  cross-section 

f„  Bluntness  factor  defined  in  Eq.  (5.4) 

D 

K  Pressure  coefficient  constant 


M 


R 

c 

SB 

S 

c 

(3 

7 

s 

A 


Mach  number 

Radius  of  spherical  cap 

Radius  of  the  base  of  a  spherical  cap  segment 
Radius  of  the  base  of  a  cone 
Planform  area  of  a  spherical  cap  segment 
Planform  area  of  a  truncated  cone 

=  n/m-'-i 

Also  angle  of  sideslip 
Gas  constant 

Angle  between  the  free  stream  floir  and  the  plane  tangent  to  the  surface 
at  a  point.  Also  compression  angle 

Expansion  angle 

Leading  edge  sweep  angle 


Subscripts 

B  Refers  to  the  bottom  half  of  a  body 
Also  the  blunted  portion  of  a  body 

C  Refers  to  the  cone  portion  of  a  body 

BC  Refers  to  a  blunted  cone  body 

I  Refers  to  the  top  half  of  a  body 
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1.0  Introduction 

Much  of  the  basic  groundwork  for  the  application  of  simple  Newtonian 
theory  to  the  determination  of  the  forces  and  moments,  and  stability  derivatives 
of  various  wing  and  body  shapes  is  presented  in  the  "First  Quarterly  Technical 
Report".  It  has  been  demonstrated  there  that  reasonably  simple  analytic  equations 
are  possible  for  all  the  important  aerodynamic  stability  derivatives  as  a  function 
of  angle  of  attack. 


In  the  "First  Quarterly  Technical  Report",  the  methods  developed  were 
applied  to  the  COne.  Only  the  static  forces  and  moments,  and  the  static 
derivatives  were  determined  during  the  first  quarter.  During  this  second  quarter, 
the  methods  have  been  applied  to  the  determination  of  the  important  damping  deri¬ 
vatives  of  the  cone,  and  the  results  are  included  in  this  report. 

Several  other  shapes  have  been  looked  at  in  the  same  detail.  Analytic 
expressions  have  been  derived  for  the  force  2nd  moment  coefficients  and  stability 
derivatives  (both  static  and  dynamic)  of  the  delta  planform  with  diamond  shape 
cross-section.  Similar  results  have  been  obtained  for  the  segment  of  a  spherical 
cap.  One  of  the  objectives  has  been  to  apply  the  spheric2LL  cap  results  to  the 
determination  of  the  stability  derivatives  of  a  blunt  cone.  Some  work  has  also 
been  done  on  the  determination  of  static  forces  2md  moments  and  their  derivatives 
for  the  delta  planform  with  elliptic  cross-section. 

Several  methods  are  outlined  for  the  determination  of  analytic  equations 
for  the  longitudinal  force  coefficients  of  blunted  surfaces.  These  methods  have 
been  applied  to  the  blunt  cone.  The  use  of  K>unt  bodies  in  the  hypersonic 
region  is  often  mandatory,  and  the  stability  derivatives  of  such  surfaces  are 
therefore  of  special  interest. 

Pressure  distributions  computed  using  simple  Newtonian  theory  are 
generally  applicable  only  for  a  limited  class  of  aerodynamic  shapes  and  limited 
Mach  number  and  surface  angles  of  attack  (m8>4.G).  A  beginning  has  been  made 
into  an  investigation  of  other  hypersonic  theories  for  computing  surface  pres¬ 
sures,  The  purpose  is  to  extend  the  methods  developed  for  computing  Newtonian 
stability  derivatives  to  lower  Mach  numbers  and  angles  of  attack. 


1 


R-63 -011-109 


Throughput  this  report  reference  will  be  made  to  equations,  sections, 
and  figures  in  the  "First  Quarterly  Technical  Report",  An  equation,  section 
or  figure  number  followed  by  the  letter  F  refers  to  the  equation.,  section,  or 
figure  in  the  "First  Quarterly  Technical  Report", 

2,0  Analytic  Equations  -  The  Circular  Cone 

2.1  General  Discussion 

The  general  methods  developed  for  determining  analytic  equations  for 
the  stability  derivatives  were  applied  to  the  cone  in  the  "First  Quarterly 
Technical  Report".  Luring  the  first  quarter,  only  the  static  coefficients  and 
static  derivatives  had  been  determined.  The  damping  derivatives  have  since  been 
evaluated  and  are  presented  in  this  report. 

At  the  outset,  it  should  be  stated  that  the  stability  derivatives  are 
computed  under  the  special  conditions  of  zero  11  (  <p  =  0)  and  zero  sideslip 

{  (3  -  0)  of  the  body.  The  analytic-  equations  rived  for  tlu  damping  derivatives 

are  therefore  a  function  of  only  angle  of  attack  for  a  gven  stipe.  For  the 

present,  the  derivatives  are  computed  assuming  simple  rytt  >n  .  theory.  prelim¬ 
inary  consideration  of  Mach  number  effects  is  discussed  in  >>ct  on  6,0. 

For  each  damping  derivative, the  damping  of  a  half  body  with  ft  at  side  up 
at  positive  angles  of  attack  is  fi'st  determined  Tnese  results  are  then  gener¬ 
alized  to  apply  to  full  bodies  at  all  angles  of  alu.ck, 

2.2  Damping  in  fitch  -  0 

ll 

As  shown  in  Section  6. IF,  the  equation  of  a  cone  in  rectangular  co¬ 
ordinates  can  be  expressed  as 

f  (x,  y,})  c  barf'd)  * *-  y*  -  *  O  (6- IF) 

The  axes  system  used,  and  the  positive  direction  of  forces  moments,  and  angular 
rotations  are  indicated  on  Figure  (1.0F) 
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As  shovm  in  Section  5. 2F,  the  damping  in  pitch  of  any  body  can  be 
expressed  in  one  of  two  ways 

Car  *  ( Cdsc<  +  {km?)*  srmn  (5.17F) 

Cm^  =  sin  (Gm^-hx)  (5.18F) 

The  constants  (K  ),  and  (K  )„  are  determined  from  Eqs.  (5.15F)  and  (5.16F) 
mq  1  mq  d 

by  integrating  over  the  body  surface.  The  partial  derivatives  appearing  in 

these  equations  are  determined  from  the  equation  of  the  surface,  Eq.  (6. IF). 

The  constants  m  and  Q  in  Eq,  (5.18F)  are  defined  by  Eqs.  (5.19F)  and  (5.20F), 
q  mq 

and  tire  repeated  here  for  convenience. 


(5.19F) 


-/ 

sin 


(5  OOF) 


The  constants  (K  ),  and  (K  ),.  have  been  determined  for  a  half 
mq  1  mq  2 

circular  cone  with  flat  side  up  under  the  assumption  that  the  angle  of  attack 
is  positive  so  that  the  complete  bottom  always  "sees"  the  flow. 


(2.1) 
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{Km%)z  =  -  t[/-  f(4)^(4f]-  #(4)]  ^ 


-  f[/--2(4)2+  -^JK4X4)^a'7  ^  sirfe*- 


-  §  tarfd  smze  (2.2) 

In  these  expressions,  ©  is  the  cone  half  angle,  and  X^  and  z^  are  the  horizontal 
and  vertical  C.G.  location  along  and  perpendicular  to  the  cone  axis  respectively. 

The  damping  in  pitch  (Cm^)  of  the  half  cone  is  based  on  the  planform  area  of  the 
cone  and  the  cone  length  (c).  The  symbols  used  here  are  the  same  as  those  in  the 
"First  Quarterly  Technical  Report",  and  one  can  refer  to  this  report  for  definitions 
of  terms. 


For  any  given  cone  angle  (  Q  )  and  C.G.  location  (x  ,  zq)  it  is  possible 

to  evaluate  (K  )..  and  (K  ).,  numerically.  The  damping  in  pitch  as  a  function 
mq  l  mq  _ 

of  angle  of  attack  can  then  be  determined  from  Eq.  (5.17F),  A  more  useful  ex¬ 
pression  for  the  damping  in  pitch  is  given  by  Eq.  (5.18F).  The  constants 

and  0  can  be  evaluated  numerically  from  (K  ) ,  and  (K  ),.,  It  thus  becomes 
mq  mq  .1  mq  2 

evident  that  as  long  as  all  of  the  conical  surface  "sees"  the  flow,  it  is  possible 
to  derive  a  simple  analytic  expression  for  the  damping  in  pitch.  It  is  interesting 
to  note  that  based  on  Newtonian  theory,  all  shapes  can  be  cast  in  the  same  analytic 
fora,  and  their  relative  damping  qualities  compared. 

It  would  be  very  useful  if  m  and  ©  in  Eq.  (o,18F)  could  be  evaluated 

q  mq 

analytically  in  terms  of  the  geometric  parameter  of  the  surface  {  Q  ),  and  the 
C.  G,  location  (X  ,  z^i.  Examination  of  Kqs.  (2.1)  and  (2.2)  substituted  in 
Fqs.  (o.ltlF)  ar.d  (S.20F)  indicates  that  the  analytic  expressions  will  be  very  in¬ 
volved  and  therefore  of  little  praetica 1  use.  Considerable  simplification  is 
possible  if  one  limits  the  discussion  to  cones  of  small  angle.  It  is  then 
necessary  to  retain  only  the  first  terms  in  the  expressions  for (K  ) .  and  (h  )., 
llamping  derivatives  for  small  cone  angles  will  be  discussed  shortly. 

It  becomes  evident  from  examining  Vqs.  (2.1)  and  12,2)  that  the  damping 
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in  pitch  is  a  rather  involved  function  of  C.G.  location.  No  simple  relationships 
can  he  established  for  transferring  the  damping  in  pitch  from  one  C.G.  to  another, 
as  is  the  case  for  the  static  derivatives  of  the  surface. 

The  pitch  damping  constants  in  Eqs.  (5.18F)  for  a  half  cone  are  plotted 
as  a  function  of  cone  half  angle  and  C.G.  location,  both  Vertical  and  horizontal, 
in  Figs.  1,  2,  and  3.  It  becomes  evident  from  the  figures  that  the  damping  in 
pitch  is  an  especially  strong  function  of  horizontal  C.  G.  The  figures  indicate 
that  the  damping  is  large  about  the  nose  and  small  about  a  point  2/3  of  the  cone 
length  aft  of  the  nose.  As  long  as  the  C.  G.  remains  in  the  cone,  the  damping 
is  not  a  strong  function  of  vertical  C.  G,  location  except  at  large  cone  angles 
and  horizontal  C.  G.'s  aft  of  the  2/3  point.  The  v.'tical  C.  G. 's  chosen 
are  rather  extreme  and  represent  the  limit  the  C.  G.  can  move  vertically  and  still 
remain  within  the  envelope  of  the  cone. 

It  is  instructive  to  look  at  the  damping  of  a  half  cone  under  the 
assumption  that  the  cone  angle  is  small,  and  the  vertical  C.G.  location  (z^/c) 
is  also  small.  This  will  be  true  if  the  C.G.  is  located  inside  the  cone  en¬ 
velope  so  that  J*  /x^j  <tan  ©.  Retaining  up  to  first  order  terms  of  small 

quantities,  Fqs.  (2,1)  and  (2.2)  reduce  to 


s  ~  ^[~4  (if)  +  (if)  ] 

(2.3) 

4  (if)  +  2  ] 

(2.-1) 

Substituting  in  Kqs,  (n, 19F)  and  (5.20F)  and  retaining  up  to  first  order  terms 
in  ©  results  in  the  following: 

<%*)*] 


,  whore 


\2.C) 
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Thus  for  a  half  cone  of  small  semi-angle  (  9  ),  the  damping  in  pitch  equation, 

Eq.  (5.18F),  is 

\  <%)*  «•■»> 

It  is  interesting  to  note  that  the  damping  in  pitch  of  a  slender  cone  is  only 
a  function  of  longitudinal  C.  G.  location.  Eqs.  (2.5)  and  (2.6)  should  be  com¬ 
pared  to  the  exact  results  of  Figs.  1,  2,  and  3.  Eq.  (2.7)  will  give  reasonably 
good  results  up  to  cone  angles  of  about  8  degrees.  It  is  possible  to  derive 

and  9  with  higher  order  terms  in  small  quantities.  This  will  increase  the 
inq 

range  of  validity  of  the  pitch  damping  equation  and  also  show  the  effects  of 
vertical  C.  G.  location. 

Fq,  (2.7)  indicates  the  strong  effect  of  longitudinal  C.  G.  location 
on  the  damping  in  pitch  of  a  cone.  Eqs.  (2.3),  (2.4),  (2.5),  and  (2.6)  are 
plotted  as  a  function  of  longitudinal  C.  G.  position  in  Fig.  4.  It  becomes 
obvious  that  the  damping  varies  widely  with  C.  G. ,  and  is  a  maximum  for  the 
C.  G.  at  the  nose  and  a  minimum  for  the  C.  G.  at  2/3  of  the  cone  length  aft  of 
the  nose. 


Vp  to  now,  the  discussion  has  been  concerned  with  the  damping  in 
pitch  of  a  half  cone  with  flat  side  up  at  positive  angles  of  attack.  The  de¬ 
rivation  of  analytic  expressions  for  the  damping  in  pitch  of  a  complete  cone  at 
all  angles  of  attack  till  now  be  discussed  in  same  detail. 


In  order  to  determine  the  damping  in  pitch  of  the  top  half  of  the 
core,  it  necessary  to  evaluate  the  constants  IK _ ),  and  (K  )s>  of  Eq.  (3.17F) 

4-  me 

for  the  top  half.  This  may  be  done  following  the  integration  procedure  outlined 
*.n  Section  o.2F  .  Ihe  method  is  the  same  one  used  for  the  bottom,  only  the 
limits  of  integration  are  different.  These  constants  will  be  valid  at  negative 
angles  of  attack  for  a  half  cone  with  flat  side  down.  Following  this  procedure, 
the  constants  for  the  top  of  the  cone  become 
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=  ~ '  I  (%■)  +  *(%*•)  }*»* 6+  x  (%‘)[/_  2  (-^P-)] s//7^  * 

*(%-)*  ]tai?0  5/^  r  =§f  (-^-)  ta.r? 9 smx0  - 

-frtan3G  sinz9  (2.8) 

(**»*)*  -  #['-  T  (%*)  +  4%f]  "  #  *<-§-)(/-  f  (#)]  + 

+  f  ['~*(^f+*(4r)*]s/^  -^^(^X#-)  **»**»*»+ 

+  §  t*n*9  sin*0  (2.9) 

These  constants  are  identical  to  those  for  the  bottom  half  of 

the  cone,  except  for  some  important  differences  in  signs  of  the  various  terms. 

Inspection  shows  that  (K  ),  for  the  top  can  be  obtained  from  (K  ),  for  the 

raq  X  niq  1 

bottom  by  substitution  -a  for  a  .  The  second  constant,  (K  ),,  of  the  top  can 

'  o  o  ’  nq  2 

be  obtained  from  (K  ).,  of  the  bottom  by  first  substituting  -a  for  a  ,  and 
mq  s.  '  o  o 

then  changing  the  sign  of  every  term. 

The  constants  m  and  0  for  the  top  are  obtained  by  substituting 

q  mq 

(K  ),  and  (K  ),,  in  Fqs.  (5.19F)  and  (5.20F) .  As  long  as  the  top  and  bottom 
wq  x  fliq  — 

of  the  cone  completely  see  the  flow,  the  damping  in  pitch  of  the  full  cone 
can  be  written  as 


=  «  e  +*x)  +  sin  ■*  <*) 


(2.10) 
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The  subscripts  B  and  T  refer  to  the  bottom  and  top  of  the  cone 
respectively. 

The  @mqg  and  fitoq^  shown  in  Eq.  (2.10)  are  not  in  general 

acute  angles  as  determined  from  Eq.  (5.20F) .  Methods  of  obtaining  the  acute 

angle  form  are  discussed  in  some  detail  in  Section  5.2F.  In  general,  for  the 

bottom  half  of  the  cone,  (K  ),  and  (K  )„  are  both  negative.  For  the  top 

mq  i  mq  i 

half  of  the  cone,  it  follows  from  the  previous  discussion  that  in  general 

(K  )-  is  negative  and  (K  )„  is  positive.  Under  these  conditions  it  be- 
mq  x  mq  i 

comes  apparent  that  ^mqg  is  "the  third  quadrant,  and  ^mq^ 
is  in  the  fourth  quadrant.  In  terms  of  acute  angles,  the  angles  for  the 
bottom  and  top  can  then  be  written  as 


9mhi  = 

*  Z7T  - 


where  the  primed  values  are  acute  angles.  By  substituting  in  Eq.  (2.10),  the 
damping  in  pitch  of  the  full  cone  becomes 


Cm*  -  -  s«»(^a *«)  -  (m$)rStn  r-  °0 


It  should  be  of  interest  to  note  that  Eqs.  (2.10)  and  (2.11)  will 
apply  whether  or  not  the  top  and  bottom  surface  are  of  the  same  shape.  It  is 
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only  necessary  that  the  constants  for  the  top  shape  and  bottom  shape  be  computed 
separately. 

If  the  C.G.  is  located  in  the  horizontal  plane  of  symmetry  of  the  cone 
(z  *  0),  then  the  damping  of  the  complete  cone  assumes  an  especially  simple  form 
since 


K)e  =  (»»$T  * 


Eq.  (2.11)  now  becomes 


-  -  +  S/h  (©/*)£  -  cx)j  (2,12) 

This  equation  could  have  been  arrived  at  from  simple  physical  reasoning  since  it 
says  that  the  pitch  damping  of  the  top  of  the  cone  at  negative  angles  of  attack 
is  no  different  than  the  damping  of  the  bottom  at  positive  tingles  of  attack  if  the 
C.G.  is  in  the  horizontal  plane  of  symmetry. 

The  first  and  second  terms  of  Eq.  (2.12)  apply  to  lower  and  upper  sur¬ 
face  respectively.  The  maximum  positive  and  negative  angle  of  attack  for  which 
this  equation  is  exactly  true  is  determined  as  the  minimum  positive  and  negative 
angles  of  attack  for  which  the  normal  component  of  free  stream  is  always  inward 
at  any  point  on  the  complete  cone  surface.  For  the  damping  in  pitch,  this  can 
be  determined  from  Kq.  (5.12F)  for  the  condition  ('\T*  f>  )  >  0.  But  this  con¬ 
dition  is  no  different  than  that  for  the  static  forces  and  moments  and  their 
derivatives  .is  determined  from  Fq.  (2.12F).  n'hen  Kq.  (3.12F)  is  applied  to  Uie 
complete  cone,  the  angle  of  attack  range  for  which  Kq.  (2.12)  is  exactly  valid 
is  -  G  <  <9  < 

It  will  now  prove  instructive  to  examine  Kq.  (2,12)  m  terms  of  a 


specific  example  and-  determine  the  relative  magnitude  of  the  two  terms  at  angles 
of  attack.  If  one  selects  as  an  example  a  5  degree  cone,  Eq.  (2.12)  is  exactly 
valid  for  the  angle  of  attack  range  -5°  £  o(  <  5°.  Also  for  a  5°  cone  it  is 
Sufficiently  accurate  to  Use  the  small  cone  angle  form  of  the  pitch  damping 
constants,  Eq.  (2.5)  and  (2.6) .  Eq.  (2.12)  now  becomes 


(2.13) 


At  plus  5°  angle  of  attack  Eq.  (2.13)  is  still  valid  for  a  5°  semi-angle  cone, 
and  the  relative  magnitude  of  the  contribution  of  the  top  and  bottom  to  the 
damping  Can  be  determined  from  the  following  ratio. 


s,n  ~  0f)  _  S/n  (s.9o°-£°\  _  sin  (.90*) 
Sin  s/n  {s.90 °+s°)  ~  s/n(/a.9o°) 


.  Q8Z5 


Thus  at  5°  angle  of  attack  the  contribution  of  the  top  of  the  cone  is  only  about 
8  percent  of  the  bottom.  At  5°  negative  angle  of  attack  just  the  reverse  is 
true,  the  bottom  is  only  8  percent  of  the  top. 


Obviously  only  a  maximum  of  8  percent  error  is  introducted  if  the 
contribution  of  the  top  is  neglected  after  5  degrees  angle  of  attack  and  the 
bottom  is  neglected  after  angles  of  attack  less  than  -5°.  The  error  will  be 
considerably  less  if  the  top  and  bottom  terms  are  retain,  d  until  the  angles  of 
attack  are  equal  to  the  absolute  value  of  the  "effective"  surface  angle  -J?  0 
For  the  5°  cone  this  effective  angle  is  5.90  degrees.  These  same  arguments  hold 
for  even  larger  cone  angles  and  a  wider  variety  of  shapes  than  the  simple  cone. 
On  the  basis  of  these  and  similar  results  it  is  possible  to  write  reasonably  ac¬ 
curate  results  for  the  damping  in  pitch  throughout  the  angle  of  attack  range 
from  Eq.  (2.12)  by  dropping  appropriate  terms,  i.'e  can  thus  write , 
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~'&rng  —  —  9  mg. 


Crr,  -  \sin{d'm^  **)  +  Sin 


(2.15) 


<*  2  0 


C/^  -  -  S//7  (&/X6  S-tx) 


(2.16) 


For  any  full  cone,  the  constants  in  Eqs.  (2.14),  (2.15),  and  (2.16)  can  be  determined 
numerically  as  discussed  previously.  For  a  slender  cone  the  analytic  forms  of  the 
constants  are  especially  simple  and  the  damping  of  a  full  slender  Cone  throughout 
the  angle  of  attack  range  can  be  written  as  follows: 

ck  <  -  —,9 


C„  =  ~  §-(/-  $-(-£•)  +  zfet] sm  if*  -  *) 


(2.17) 


3jfo  <  cx  <. 


CmQ-  -  y['-  I  (yr)  +  *(!?-)  }[*”(¥*  +«)  +  s">  (if*-*)}  (2.18) 


«13-Z-9 


Cm  [/-#  (£■)  ♦  4%}*]  *«  (?***) 


(2.19) 


The  damping  in  pitch  at  zero  angle  of  attack  for  a  slender  cone  becomes 
from  Eq.  (2,18), 


(2.20) 


If  the  C.G.  is  not  located  in  the  horizontal  plane  of  symmetry  (*  /  0) , 


U 
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then  the  constants  for  the  top  and  bottom  of  the  cone  are  not  identical  and  the 
more  general  form  of  the  damping  in  pitch  equation  of  the  Complete  cone,,  Eq.  (2.11) , 
must  be  used.  The  constants  for  the  top  and  bottom  must  then  be  evaluated  separate¬ 
ly.  From  Eq.  (2.11) ,  good  approximate  expressions  for  the  damping  in  pitch  through¬ 
out  the  angle  of  attack  range  can  be  derived  using  the  same  procedure  that  was 
used  for  the  slender  cone.  From  Eq.  (2.11),  these  equations  become 


<X  <L  —  & 

Crr>%  =  -  ~  «)  ( 2 • 2i) 

**  —  &m$-r 

C m  -  ~  (m$)e  s,n(&f*yB+  •*)■  (m?>T  s,n  )  (2.22) 


Z  r 

C^=-(^)e  (2*23) 

Although  the  discussion  here  has  been  concerned  with  the  full  cone,  the 
form  of  the  results  will  be  the  same  for  any  shape  surface  in  Newtonian  flow.  In 
fact,  the  upper  and  lower  surface  need  not  be  the  same  shape  for  Eqs.  (2.21),  (2.22), 
and  (2.23)  to  be  applicable.  It  is  only  necessary  that  the  constants  in  these 
equations  be  determined  separately  for  the  top  and  bottom  shape.  Similar  results 
for  the  damping  in  pitch  have  been  determined  for  other  shapes  and  these  will  be 
presented  in  subsequent  sections. 

The  other  damping  derivatives  can  be  determined  following  essentially 
the  same  procedure,  and  these  will  be  discussed  in  the  sections  that  follow. 


V 
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2.3  Damping  in  Yaw  -  Cnr 

It  is  demonstrated  in  Section  5.3F  that  the  damping  in  yaw  of  an  aero- 
dynamic  shape  can  be  expressed  as 


Cffr  -  ( kn* ),  COS*  +  {Knriz  Sjn<x  (5. 28F) 

Cnh-  nr>nr  sin  {&nr  r«)  (5.31F) 


where 


(5.32F) 

(5.33F) 


The  constants  (  K  and  (  K  )_  are  determined  from  Eqs.  (5.29F),  and 

nr  nr  c. 

(5.30F)  using  the  equation  of  the  surface,  Eq.  (6.1F) . 


For  a  half  cone  with  flat  side  up  at  positive  angles  of  attack,  the 
foil owing  expressions  for  the  constants  have  been  evaluated. 


[KnrU  -  -  IT  [  /-  4  +  tAn  Q  - 

+  s.;nz&  ~tr  tas?&  (2  2H 

i  [/-  i&h 

For  any  com'  half  angle  ( 8  )  and  C.6.  location  (x^,  ,v  ),  it  is  possible  to 
evaluate  numerically  (  X  ),  and  {  X  and  also  the  constants  *  and 

jo  x .Fqs.  (5.23F)  and  to.diF)  will  then  give  the  damping  in  yaw  as  a  function  oi 
angle  of  attack  for  positive  angles  of  attack. 

Id 


R-63 -011-109 

The  damping,  in  yaw  is  a  function  of  both  longitudinal  and  lateral  C.G. 
location.  No  simple  C.G.  transfer  relationships  can  be  derived  as  they  were 
for  the  static  derivatives. 

The  yaw  damping  constants  in  Eq.  (5.31F)  ,with  the  acute  angle  form  for 
,  are  plotted  for  a  half  cone  as  a  function  of  cone  angle  and  C.G.  location 
in  Figs.  5 ,  6,  and  7.  The  yaw  damping  is  an  especially  strong  function  of  longi¬ 
tudinal  C.  G.  position  and  varies  little  With  lateral  C.  G.  location  for  moderate 
Cone  angles,  &  <15°.  The  damping  is  large  about  the  nose  and  small  about  a 
point  2/3  of  the  cone  length  from  the  nose.  Only  lateral  C .  G.  locations  that 
lie  within  the  Cone  envelope  have  been  Considered  in  these  figures,  |yo/*0|  <- 
tan  Q  . 

It  is  interesting  to  look  at  the  damping  in  yaw  under  the  assumption 
that  the  cone  half  angle  is  small,  and  the  lateral  C.  G.  location  is  small , 
|yQ/xo|  <  tan  O  .  Retaining  only  up  to  first  order  terms  in  small  quantities, 


the  constants  can  be  written  as 

[Knr),  =-»•[/- 4  (2.26) 

(Me  =  -  l[' -T  (if)  *  <2-2T> 

mor  -  +  (2-28) 

Qnr  z  ir  +  e'nr  >  ''hero  Qnr  =.  (2.29) 

For  sufficiently  small  cone  angles,  Eq.  (d.dlF)  thus  becomes 

(& &  +«)  (2.80) 


Kqs.  (2.28)  and  (2,29)  should  be  compared  to  the  exact  solutions  of  Figs,  5,  6, 
and  T,  These  comparisons  indicate  that  Eq.  (2,30)  will  give  reasonably  good 

11 
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results  for  Gone  semi-angles  less  than  about  10  degrees.  The  strong  effect  of 
longitudinal  G.  G.  location  on  the  damping  in  yaw  io  indicated  by  the  plot  of 
EqS .  (2.26)  through  (2. 29)  in  Fig.  8.  The  figure  indicates  that  the  damping  in 
yaw  is  a  minimum  for  the  C.G.  located  at  2/3  Of  the  Cone  length  from  the  nose. 

The  damping  in  yaw  of  the  full  Cone  is  now  fairly  easy  to  evaluate  for 

all  angles  Of  attack.  If  the  constants  (K  )  and  (  K  )  -  for  the  top  half 

nr  l  nr  2 

of  the  full  cone  are  determined ,  it  will  be  found  that  the  (  K  )„  of  the  top 

nr  1 

and  bottom  half  are  the  same*  and  the  (  Knr  of  the  top  is  opposite  in  sign  to 

the  (  Kjjj.  ).  of  the  bottom.  This  same  relationship  existed  for  the  (  K  )., 
c.  mq  J- 

and  (  K  )g  of  the  top  and  bottom  with  the  C.G.  in  the  horizontal  plane  of  sym¬ 
metry.  It  therefore  follows  that  the  damping  in  yaw  for  the  complete  cone  is 
analogous  to  Eq.  (2.12)  and  can  be  written  as 


-  m 


nr 


[s/n  {&'nr-  +«)  +  sin  (&'nr  -<*)j 


(2.31) 


The  damping  in  yaw  throughout  the  angle  of  attack  range  can  now  be  determined 
by  simple  analogy  to  the  damping  in  pitch  equations. 


c*  —  ~  3  nr 


Cnr-~^nr  S'r'ie'nr--*) 


(2,32) 


Cnr~~^>nr  [*">  (<?hr  +  «)  +  *">(*«»-  *  *)]  (2.33) 

<*  ^ 


Cn^--'»nr 


(2.3-1) 


U' 
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If  the  upper  and  lower  surface  are  not  of  the  same  shape,  then  the 
constants  for  the  upper  and  lower  surface  must  be  determined  separately  as 
indicated  in  Section  2.2  when  discussing  the  damping  in  pitch.  In  such  a  case, 
the  yaw  damping  equations  will  be  analogous  to  the  pitch  damping  Eqs.  (2.21) , 
(2.22) ,  and  (2. 23) ,  and  they  can  be  written  as  follows: 


64  ~  ~&nre 


Oir  ~  ~  ^nrr  s,n  "  ®0 


(2.35) 


-e'nrB  e'r>rT 

Cnr  =  ~  +  «)  “  "Wr  s>" 


(2.36) 


9/>rr 


Cn^~mnre  S">(*»ra+'X) 


(2.37) 


The  subscripts  T  and  B  refer  to  the  constants  for  the  top  and  bottom  suiface 
shape  respectively. 


For  a  full  slender  cone  Eqs.  (2.32),  (2.33),  and  (2.3*l)  assume  an 
especially  simple  form. 


04  *  -f6 


c„  =  - 1  ['- 1  +■  j?  - «) 


(2.38) 


H' 
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37T 

-f 


&  <  cx  < 


sir 


6 


[*»{£:*«)*  <2.39) 


<X 


c„*  -  f  [/- 1  (■%•>  -  *&f 


(2.40) 


At  zero  angle  of  attack,  the  damping  in  yaw  of  a  slender  cone  can  be  derived  from 
Eq.  (2.39) .  The  result  is  identical  to  the  damping  in  pitch  at  zero  angle  of 
attack. 


(2.41) 


2.4  Damping  in  Roll  -  C. 

XP 

From  Section  5. IF 


flow  can  be  expressed  as 


the  damping  in  roll  of  a  surface  assuming  Newtonian 


where 


/* 


(5.43V) 


Cjt  =■  *"*(&*/>  f«) 


(5.44F) 


(5.45V) 


I 


i'V*i 


<5.4t?V) 


.IT 


-  t*n 
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Using  the  equation  of  the  surface,  the  constants  and  (K^p)2  can  be 

determined  as  indicated  in  Section  S.4P.  For  a  half  cone  with  flat  side  up, 
at  positive  angles  of  attack,  these  constants  become 


*  (%f  ]  s,n  &  tos  & 


(2.42) 


'**■ 


cos2 o 


(2.43) 


For  any  cone  angle  ( Q  )  and  C.G.  location  (y^,  zq) >  the  constants 
can  be  evaluated  numerically,  and  the  damping  in  roll  of  a  half  cone  can  be 
determined  as  a  function  of  angle  of  attack  from  Eq.  (5.44F).  It  is  interest¬ 
ing  to  note,  that  with  the  C.G.  on  the  cone  axis  (yQ=  zq  =  0),  the  damping  in 
roll  is  always  zero,  as  one  would  expect. 

For  those  cases  where  the  lateral  displacement  of  C.G.  is  zero, 
yo/c  =  0,  the  constants  assume  an  especially-  simple  form 


(Kjf  )t  -  ~*w  s'°  #  c*s  & 

~  ttr  &  e  1 tan  G  S tn*&  (2.44) 

(*> A  = 

=  -  f  (-^-f  Sm*G  (2.45) 

If  the  further  assumption  is  made  that  the  cone  angle  is  small,  and  only  up  to 
first  order  terms  in  $  are  retained,  then  the  constants  assume  the  following 
values ; 

18 
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U&b  =  0  (2.46) 


(2.47) 


§_ 

3 


(2.48) 


he  *  7r+  '  where  (2.49) 

The  roll  damping  constants  determined  using  Eqs.  (2.44)  and  (2.45) 
are  plotted  as  a  function  of  cone  angle  in  Fig.  9.  These  results  compare  quite 
well  for  &  <  10°  to  the  small  angle  expressions,  Eqs-  (2.48)  and  (2.49). 

For  a  slender  half  cone,  the  damping  in  roll  as  a  function  of  angle  of 
attack  and  vertical  C.G.  position  thus  becomes 


Cf  .  *  -  I (^f  sv'i***  +*)  (2.50) 

/ 

The  damping  in  roll  of  the  full  cone  is  easy  to  determine  once  the 
constants,  (K^p)^  and  for  the  too  half  of  the  cone  are  evaluated. 

Analysis  shows  that  these  constants  con  be  derived  from  those  for  the  bottom 
half  by  first  substituting  -  for  c  and  then  changing  the  sign  of  all  the 
terms  in  (K  )*,  .  It  will  be  remembered  that  the  same  procedure  applied  to 

the  damping  in  pitch  constants.  By  analogy  to  Eqs. 12.21),  (2.22),  and  (2.23), 
the  damping  in  roll  equations  become 


•  a 
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’f  ~j*fe 


CAr*'m-tfr^k"x) 


(2.51) 


Vtfg-*-  6Jfr 


Cer  *  -»jfs  ^fe  +«-)-"Ufrs'* 


(2.52) 


(2.53) 


The  subscripts  T  and  B  refer  again  to  the  constants  for  the  top  and 
bottom  surface,  and  these  constants  are  generally  different  because  of  lack  of 
symmetry  of  the  upper  and  lower  surface  with  respect  to  a  C.G.  outside  the  hori¬ 
zontal  plane  of  symmetry.  Of  course  these  same  equations  will  apply  if  the  upper 
and  lower  part  of  the  body  are  different  in  shape  provided  the  constants  are 
evaluated  for  each  shape.  In  the  case  of  the  cone,  the  damping  in  roll  constants 
happen  to  be  an  even  function  of  the  vertical  and  lateral  position  of  the  C.G., 
see  Eqs.  (2.42)  and  (2.43),  and  therefore  the  constants  will  be  identical  for  the 
top  and  bottom. 

For  the  slender  cone,  using  the  results  of  Eqs.  (2.48)  and  (2.49),  the 
damping  in  roll  of  the  full  cone  assumes  an  especially  simple  analytic  form 
throughout  the  angle  of  attack  range. 


20 
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c*  <  -  Sfg 


(2.54) 


-2Z&  <  (X  <  ^  9 


3ir . 


Cy-  -  f(W  M *fe+«) *  *"> HP* «*)  ] 


(2.55) 


of  >  ^9 


Cj^  -  -  §-($r)  S/n  9  +  <x) 


(2.56) 


The  damping  in  roll  for  a  slender  cone  at  zero  angle  of  attack,  with 
no  lateral  C.G.  displacement,  can  be  obtained  from  Eq.  (2.55) 


C.  -  -  4*r(±'  )\.>  (2.57) 

'/ 

2.5  damp  in;;  in  Yaw  dee  to  Rolling  -  C 

*"*  "  np 

l'he  general  analytic  form  of  the  damping  in  yaw  due  to  rolling  is 
presented  in  Section  5.5F  and  repeated  here  for  convenience. 


A 


(*,.  ), 


\'S  > 


l1 


r/f 
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cnp  -  ™nf>  M 

where 

mr>f  ~  *  {K/7/s)z 


(5.55F) 


(5.56F) 


tan 


l^Ai p)/ 

{Knfh 


-/ 

=  s/n 


_ fan-fil/ _ 

{Knf)z 


(5.57F) 


The  constants  have  been  determined  for  a  half  cone  with  flat  side 
up  at  positive  angles  of  attack  as  explained  in  Section  5.5F.  The  results  are 

(*/?/>)/  =  4  (■£■)]*"» &  +  & (-%-f  s//7^ - 

tan  0  smz0  (2.58) 

{KnfV = -  -f-  (4*-)  [/-  4  (-^)1 + **  (4)^"  <9  - 

"f  (4)(^ 3//)^  '  ^(4f # afl  ^  •S/n^  (2.59) 


For  any  cone  angle  (  Q  )  and  C.G.  location  (xq,  yQ,  zo),  the  constants 
can  be  evaluated  numerically.  The  Cn  as  a  function  of  angle  of  attack  will  then 
be  defined  by  Eq.  (5.54F)  or  (5.55F).P 

For  a  slender  cone,  assuming  y  / c  -  0,  and  retaining  only  up  to  first 
order  terras  in  the  cone  angle,  considerable  simplification  is  possible. 


•JO 
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(Knf),  =  -  *?(¥>[>-  <2*60> 


^nf)z  z  ~  $f)  [/~  I  (^)J  (2.61) 

”>nf>  =  #(#■)[/-  iity]  (2.62) 


Gnp-IT*  0hf>  where  ^  r  *2*^»  (2.63) 

Thus  for  slender  half  cone  at  positive  angle  of  attack,  Eq.  (5.55F) 

becomes 


Cnf>  ~  ~  T  (4H(/“  i  C"^)]  +  «0 


(2.64) 


The  constants  for  Cnp  have  been  computed  as  a  function  of  cone  angle 
(  Q  )t  for  fixed  longitudinal  C.G.  locations  (x0/c),  and  zero  lateral  C.G.  dis¬ 
placements  (y0/c  =  0).  The  constants  «np  and  &np  are  plotted  as  Figs.  10,  11, 
and  12.  The  results  show  that  Cn^  is  a  linear  function  of  vertical  C.G.  position, 
and  varies  widely  with  longitudinal  C.G,  position.  The  damping  is  a  minimum  for 
the  longitudinal  C.G.  located  2/3  of  the  cone  length  from  the  nose,  and  it  will 
even  change  sign  depending  on  the  C.G.  location.  The  results  also  demonstrate 
that  Eq.  (2.64)  will  give  good  results  for  half  cone  angles  less  than  about  10 
degrees. 


The  CPp  of  the  complete  cone  can  be  found  once  the  constants  (Kjjp)^ 

and  (Knp)2  °f  the  upper  half  are  determined.  These  were  evaluated  following 
the  procedure  in  Section  5.5F.  It  was  established  that  the  (K^)^  of  "the  top 
can  be  obtained  from  (K np)^  of  the  bottom  by  first  substituting  -zQ  for  zn, 
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and  then  changing,  the  sign  of  all  the  terms.  The  (Kn  of  the  top  is  obtained 

from  the  bottom  value  by  simply  substituting  -a  for  zq.  The  procedure  to  be 
followed  in  obtaining  the  Cnp  throughout  the  angle  of  attack  range  is  similar  to 
that  used  for  the  other  derivatives. 


2.6  Damping  in  Roll  Due  to  Yawing  -  r 

The  damping  in  roll  due  to  yawing  (Cg  )  for  the  Cone  is  the  same  as 
the  C-  since  the  Constants  are  identical. 

(*/y>)/  ~  (*£?)/ 

{Knp)z  ~ 

It  is  also  true  that  (K  ).  and  (K  .  )„  for  the  top  half  of  the  cone  bear 

ir  1  ir  2 

the  same  relationships  to  the  constants  for  the  bottom  half  as  that  shown  for 

(K  )-  and  (K  )„.  The  alternate  form  of  the  C»  constants,  (m  )  and  (&'.  ), 
np  l  np  i  *r  JLv  It  1 

are  also  presented  in  Figs.  10,  11,  and  12. 

3.0  Analytic  Equations  -  Delta  Planform  with  Diamond  Shape  Cross-Section 
3.1  General  Discussion 

Another  aerodynamic  shape  that  has  been  studied  in  some  detail  is  the 
delta  planform  with  diamond  shaped  cross-section  in  planes  perpendicular  to  the 
longitudinal  axis.  All  the  surfaces  of  such  a  shape  are  flat  planes.  The  im¬ 
portant  stability  derivatives,  both  static  and  dynamic,  have  been  determined  for 
this  shape.  In  the  next  quarter ,  the  shape  will  be  studied  in  some  detail  with 
controls  added  for  purposes  of  trim  at  angle  of  attack.  Stability  derivatives 
will  then  be  determined  under  trim  condition,  a  primary  goal  in  this  whole  study. 

Results  are  presented  for  only  a  half  body  with  the  flat  side  up 
at  positive  angles  of  attack  so  that  only  the  bottom  surfaces  "see"  the  flow. 

No  discussion  of  analytic  results  are  included  here  for  the  complete  surface, 
top  and  bottom,  at  all  angles  of  attack  as  was  done  for  the  cone.  The  procedure 
to  be  followed  to  obtain  the  derivatives  of  the  whole  body  is  essentially  that 
discussed  for  the  cone,  and  the  cone  results  may  be  referred  to  for  these 
details . 
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In  the  derivation  of  all  the  stability  derivatives,  the  partial  deri¬ 
vatives  of  the  equation  or  equations  of  the  surface  are  required.  The  body 
orientation  and  general  shape,,,  as  well  as  the  surface  equations,  will  be  shown 
here  for  convenient  reference. 


Surface  (1) 


=  >  +(*«/> -A-)y  +  ( cote )%  z  o 


(3.1) 


Surface  (2) 


F[xtf,^)zy.  ~(t±nX)f+  (cot  0)^-0 


(3.2) 


3,2  Longitudinal  Force  and  Moment  Equations 

The  forms  of  the  longitudinal  force  and  moment  coefficient  equations  as 
derived  in  Section  3. OF  are  repeated  below.  The  analytic  equation  for  each  co¬ 
efficient  can  be  written  in  either  of  the  two  basic  forms. 
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CA  i  \Ka)/:  cos%  +  {kA)2  smc*  +  {kA)3  s/tm  cas<x  (3.185") 

Q  =  +  *i  S"1*^**)  (3.19F) 

C*,  =  (<y)y  C05  <*  +  5//7  <*  +  (^)5  j//>7c>(  C«o(  (3.31F) 

C„-  n0  +  n,  smZ(Bn+<x)  (3 . 35F) 

-  *  +  (A'/rtie  S/n*o(  +  (a^)3  5//70f  c«jr o<  (3 . 40F) 

=  mQ  +■  m,  S//)Z(0m  +  <*)  ( 3 . 44F) 


The  integral  relationships  and  equations  for  evaluating  these  con¬ 
stants  are  presented  in  Section  3. OF.  These  integral  relationships  are  used 
in  conjunction  with  the  equations  of  the  surface,  Eqs.  (3.1)  and  (3.2). 

For  the  half  body  at  positive  angles  of  attack,  the  constants  have 
been  evaluated  and  are  shown  below. 


(**),  =  **£*£*13-2- 

/  *■  S/n  0  to.*  A 


(v  \  -  zs/n9 


(^)j  *  77 


4 


stnM0  tci*. A 


0.c~0 


a.  - 

/<* s,n*9to.»\A. 


& 
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(*!v)/  = 


t  /  sm  Q  £&n*JU 


n0  -  o 


a Zees.2' 6 


~  /+s/nz&  it**,  *--A- 


y  \  _  4  S//10  CoS  0 

N'*  ~  /  +  Srn*0  6».*,  K 


n,-- 


/+s/nx0 

9n'-e 


It/  \  -  ~  X  6  ton'd 

\Kth)l  -  /  +  *,»*■  a 


z  &  £a.n  A* 


noQ  =  o 


>  1. 
S/r>  & 


-Jrif'ft.H  9 


/r>‘  /*  s/n  *0  ttur,  ll 


/  -  ^rtCLnO  S/*7a& 

Wi  -  — JL 


li  **  /  +sin*0£tLnx-A~ 


0fnz%  +  #m  > where  ^  -  <$* 


The  more  useful  form  of  the  analytic  expressions  for  the  longitudinal 
force  and  moment  coefficients  are  expressed  by  Eqs.  (3.19F),  (3.35F),  and 
(3.44F).  In  ter  j  of  the  evaluated  coefficients  these  equations  become, 


_  Z  t*.n  9  _  */  _  ,  \ 

Ca  '  - S - 3 —  sm  (  &  +  oU 

*  /+sin*9tM*JL  1 


(3.3) 


CM  -  - 4 — — -r—  Stf7*’(Q  +(*) 


(3.4) 


_  —  y-ts.^6  */  x 

C„,  =  - (0+o<) 


(3.5) 


Eq.  (3.5)  is  evaluated  for  a  C.  G.  located  2/3  of  the  body  length  aft  of  the 
nose  on  the  longitudinal  x  axis.. 

It  is  interesting  to  note  that  in  all  cases  the  "effective"  body 
angles  (  0&,  @n,  &\m)  are  equal  to  the  geometric  angle  9  .  When  these  expressions 
are  Used  in  determining  the  forces  and  moments  of  the  full  body  throughout  the 
angle  of  attack  range,  the  resulting  equations  are  exact  since  complete  masking 
or  unmasking  of  the  bottom  and  top  occurs  precisely  at  angles  of  attack  -  9  and 
4  Q  respectively.  The  complete  surface  results  can  be  obtained  in  a  manner 
identical  to  that  Shown  for  the  full  cone  in  Section  6. IF. 


3.3  Lateral  Directional  Static  Stability  Derivatives 

The  static  lateral-directional  stability  derivatives  for  a  surface 
can  be  defined  analytically  in  the  following  forms  as  demonstrated  in  Section 
4.3F. 


COSo<  +  S//7e* 

(4.28F) 

C Yp~  Y(3  s,n  (&Y/3  +e*) 

(4.32F) 

%  =  coS«  +  (***)*  3,"°< 

(4.35F) 

(4.38F) 

C^=  (A^),  cor  +  {Knfl)z 

(4.43F) 

cnp*  s;*{enf3+<x) 

(4.40F) 

The  integral  relationships  and  equations  for  evaluating  these  con¬ 
stants  are  shown  in  Section  4.3F.  The  equations  of  the  surface  that  must  be 
used  are  Eqs.  (3.1)  and  (3.2).  For  the  half  body,  the  following  values  of 
the  constants  were  obtained. 
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I  v  \  _  -4  ia.n.6.  S/HZS  ta-n  Z-A~ 

W"  /  *-s/nz  6  ieLn^^t 


y  _  4£m&  §*aJLj&uds. 


/*S/n*9  t^nZA 


(xYfl)zz 


-4.  S/r)  &jte.r,ZJL 


Oypi7f+  9^  ^  where  <2^g  =  0 


(**/?)/“ 


&  tajt  &  tan^A.  [ta-n  0  4  3  (4j£j]  ^ 


{Kfi 


/  r  s/n*ff  t0j,*LA. 

-  smQcoSdlj-  ta.n  Oian*A  [tan  04  3(%-)]^ 

/f-Stn^O  ta.n*JL 

Is"5*  {' 

t~  tan  O  t*nr_A  [tan  04  3(-%^j  ^ 

/nStA*-0ttu»*JL 

s  ”  *  *j/s 

>  where  ~  & 

^  tun  &  sin*8  z 4a.n  tt  £ f -  j 

/  ■*■  -S//7  ^  £jL.4n  jl 


{Xnfl}*  - 


^  srrf&^+tta-ntA  [/-  4  }• 


/  r  tayi 
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Cfn^  -  s*#  *  <*)  (5.18F) 

The  method  for  evaluating  the  constants  in  these  equations  is  also  shown  in 
Section  5.2F. 

For  the  half  diamond  shape  body  at  positive  angles  of  attack,  the 
constants  assume  the  following  form: 


/+  S//?*0  ta.n*-A. 


8  [(%A  y  g?)  -  (•%*)*]  tern  6  sm  *9  _ 

/  *■  St/>  ta.n*-A. 

/  +  T/o  z0  / *  stn*9 


✓  *  s/*  V  U. 


/  ^  s»/i  *<9  #*>**-A 


8  [&)*♦  ?  (-^)  -{£)*]  S»"*Q  + 

/  +  Stf»  *# 


Jffa)  [/-3 (^)] tan 6S/r)*0  -f  *"?*< 

+  /  /* 


■Jt 
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+  -j  tan  &  S/r)*&y 


&/T)q.  -If  +  i  where  &'  =  6 

The  pitch  damping  constants  are  rather  involved  functions  of  the 
longitudinal  and  verticalC.  G.  location.  These  constants  can  of  course  be 
evaluated  numerically  for  any  C.  G.  location  and  values  of  the  angles  9 
and  -A.  .  In  the  case  of  small  0's,  if  tan  JL  is  of  the  order  of  one  or  less, 

I  *0  I 

and  the  vertical  C.  G.  is  located  within  the  body  (  —  <  tan  9  ),  considerable 

i  x0  i 

simplification  is  possible.  Retaining  only  up  to  first  order  terms  in  B  and 
zo/c»  E<1*  (5.18F)  becomes 


Equation  (3.9)  should  be  compared  to  Eq.  (2.7)  for  the  half  cone.  In  the  limit, 
as  B  approaches  zero,  liq.  (3.9)  becomes  the  damping  of  a  flat  plat-  delta 
surface  at  angle  of  attack. 

3.5  Lateral-Directional  Damping  Derivatives  -  CUr,  Cc^,  Cn  ,  Cjr 

The  form  of  the  lateral -direct iona.1  damping  derivatives  are  determined 
in  Sections  5.3F,  5.  IF,  o.5F,  and  5.C5F.  The  integral  relationships  and  equations 
for  evaluating  the  constants  in  the  damping  derivative  equations  are  also  pre¬ 
sented  in  these  sections.  For  convenience,  the  two  forms  of  each  damping  de¬ 
rivative  equation  are  shown  again  below. 
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Cnr  *  (*/*-)/  ^yo<  +  (*„,-)«  s>/>* 

(5.28F) 

C/,^  =  ftl fit  S">  0nr+<x) 

(5.31F) 

CJlf)  -  faf)*  Cos*  +  ('%>*)* 

(5i43F) 

fa  -  to**  s,n  fa-#  +*) 

(5.44F) 

Cfi  -  {Krtpjt  c°s  «C  +  S/a #c 

(5.54F) 

=  mnp  s/fi  {&fip  +<*) 

(5.55F) 

fa  Cos *  +  (fay*5*** 

(5.64F) 

fa  *  "Ur  **  far  *  «) 

(5.65F) 

The  constants  in  these  equations  for  the  half  body  at  positive  angles 
of  attack  have  been  evaluated.  It  was  discovered  early  in  the  integration  pro¬ 
cess  for  determining  the  constants  that  some  simplification  was  possible.  For 
the  shape  under  consideration,  the  following  relations  hold: 


to nr  ~  ~  iec 

toXf,  =  -  sec 

ft) ftp  ~  5ec  ^ 

nUr  -  ^ far)# 


&r>r  ~  ^*r  >  ^ 

where  ®nr  ~  ® 

where  &  /p  *  6> 

%  --  4  -  *  * 
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An  examination  of  the  damping  in  pitch,  and  the  lateral-directional 
static  stability  derivatives  discussed  in  Sections  3.4  and  3.3  shows  that 
similar  relationships  hold  for  these  derivatives.  Listed  below  are  only  the 
values  Of  (KjU.)2>  (K^p)2>  ^Knp^2’  s*nce  ^roin  these  the  Other  con¬ 

stants  Can  be  easily  determined. 

[Knr)z  s  'h  ~  ^  (^)  +z  Hr~f  1  -f[/~  %  (%~) *£ (~^)  ]  u  *  T  cot-A-} 

✓  ■*»  .S//?  *0  to-n  *A. 

{lcoTJL+ ]  “  ♦(£)<*«*  - 

-  [/^^(^-)* +  cotji  +  sm*&  + 

+  *(%-)[/  +  tas>A^  ta-n  G  s/n*&+  t**JL  ttjn'O  snrfo} 

{X*f>)x.=(Kjlr)t-  /I  t-  |-Cot_X  + 

-  '$(’£•)['-*('*  i-  -^-) tan't]  - 

-  ^  [/*^(-^)  -  (-^-)(/-  ^  ^  co/ stn*&y 

It  is  of  interest  to  note  that  the  cross-damping  derivatives,  C„p 
and  C |  ,  will  he  identical  since  (X  L  and  (K  V,  are  equal. 

411-1  1C  4.  I  tv 

Considerable  simplification  is  possible  in  the  equations  for  evaluating 
the  constants  if  one  assume  that  the  lies  on  the  x  axis  <*  'c  0,  y  /c  *  t>), 

i)  X> 


{Knflz  '  -  77 s#«4«w»5t  +  T  c®^]  *"» «  - 

~  j £/-  ('^X/~  ^  fa*»-A  -*-  -j  £«*■ l/l J  ta/j  0  sin*0  y 

If  the  further  assumption  is  made  that  d  is  a  small  angle,  and  tan  JL  is  of 
the  order  of  1  or  less,  then  further  simplification  is  possible  in  the  constants. 
Retaining  only  the  lowest  order  term  in  0  for  each  constant  ve  have 

OWU  *-♦<['-*($>]»  ['-  §(*)«&)*  ]<«i*  i 

~  T  Cdf  -A 

(**?)*  ~  (**r)v?  *  *  f  ('!?)  t  -j-  ^t]  £ 


Rith  these  assumptions,  the  constants  m  ,  b,  ,  a  .  and  m.  are 

nr  fp  np  ir 

identical  to  ^|p’ *'  tK^V,,  and  <K^VA  respectively. 
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Now,  for  the  half  diamond  shape  body  at  positive  angles  of  attack,  the 
lateral-directional  damping  derivatives  assume  the  following  simplified  form1. 

%=  -*{[/-  *«&**)  (3.10) 

Cs  r-4  cot*A.  sin(e+<x)  (3.n) 

*v 

Cn  =  CJ  =^[/-f  ®  +  £catjL]0s,n(e+*)  (3.12) 

1°  t 

In  the  limiting  case  of  0  =  0,  Eqs.  (3.10),  (3.1l),  and  (3.12)  become 
the  damping  derivatives  of  a  flat  plate  delta  surface  with  sweepback  -A  .  For 
this  case 


C  -  -  %  COt*J_  S //?<*  (3.13) 

jf  ■* 


(3.14) 


It  is  possible  to  look  at  other  limiting  forms  of  the  surface  with 
diamond  shape  cross-section.  If  one  assumes  that  cot  A  tan  &  ,  then  the 

shape  is  one  that  can  be  inscribed  in  a  cone.  One  can  then  make  interesting 
comparisons  between  the  derivatives  of  this  shape  and  those  of  a  cone. 

4,0  Analytic  Equations  -  The  Elliptic  Cone 

4,1  General  Pis  ussion 

An  aerodynamic  shape  of  some  interest,  upon  which  some  initial  investi¬ 
gations  have  been  made  is  the  elliptic  cone.  Only  analytic,  equations  for  axial 
and  normal  force  coefficients  have  been  determined  up  to  now  and  the.  results  are 
presented  here,  I'he  moments  and  the  lateral-directional  stability  derivatives 


An 


will  be  determined  in  the  next  quarter. 
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The  results  for  only  a  half  body  with  flat  side  up  at  positive  angles 
of  attack  are  presented.  The  equations  for  the  complete  Surface  at  all  angles 
of  attack  can  be  easily  obtained  using  the  procedure  discussed  for  the  Cone  in 
Section  6. IF. 

A  sketch  of  the  half  elliptic  cone  body  and  the  equation  of  the  surface 
required  in  the  analysis  is  shown  below. 


Equation  of  the  surface: 


1.2  Longitudinal  Force  Equations 

The  general  form  of  the  longitudinal  force  equations  are  derived 


in  Section  3f)F  and  presented  in  Section  3.2,  The  integral  relationships  and1 
equations  required  to  determine  the  constants  in  these  equations  are  shown 
in  Section  3. OF.  For  the  half  elliptic  cone  body  the  Constants  become 


-  ffcosj.  tan  Qsin  Q 


)„  v  „  ar.CdsJL. sin  0cos9 
'  A  z~  cos-yl+  swS 


(*4)3 ~ 


/I  Z 

2  s//)  O  » 

/+  1 

]* 

s/n  0 

V  V  cos  JL  > 

3k  ^ 

/-I 

[/-tSSf-J 

f 

cosJL 

</ 


(*#)/- 


cotJL 

t*-n  0 


{*)))/  ~ 


sm  Q 


l  +  sm 

,-s.neUnAi*ggt-lt 


cot-1 

t&n& 


>/ 
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f  Km) 2  - 


-  Zcos" 9 


sm  dita-n  Ji-CoV 


/- 


Sme8  -Sec2.i 


sm& 


i^A- 


cot*6 


ten  {sin&^ta.n*A  -  cot*g  ) 


.GOtziln..  <  / 
tan& 


(^/V^  - 


2 


[" 


S//7  0[cot£0  -  tcLn*jC)  I  2.  s/ri  Q  V cot *0  -  ta.n  *iA. 


It  s/n8  ^jcoity  ~  t^JL 

/  -  S/n0  )jcot*0-  tasi^A. 


CotJL  ~  , 
ter,  6 


{^1)3  ~ 


Z.fr  cos  & 
CSC  9  *  sec -A. 


For  any  values  of  the  angles  .A  and  Q  ,  numerical  evaluation  of 
the  constants  in  Eqs.  (3.18F),  ( 3 , 19F ) ,  (3.31F),  and  (3.35F)  become  obvious 
following  the  procedure  outlined  in  Section  3. OF. 

It  should  be  of  interest  to  examine  slender  elliptic  cones  of 
various  fixed  ratios  of  the  two  axes  of  the  ellipse  From  the  previous  sketch, 
the  ratio  of  horizontal  to  vertical  axes  cf  the  ellipse  ( e  )  is  obviously 

e  - 

tan  8 

If  the  assumption  is  now  made  that  the  angle  9  is  small,  and  e  is  of  the 
order  of  one,  it  is  possible  to  obtain  quite  simple  analytic  expi'essions  for 
the  constants  in  terms  of  e  and  the  angle  6  .  If  the  value  of  e  is  fixed, 
then  the  constants  become  only  a  function  of  &  ,  For  small  <9 ’s,  the  con¬ 
stants  are  listed  below  for  three  ellipses. 
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e  =  2,0 


a 

1.8642 0 

0 

ai 

2.0944  6* 

K 

1.452  9 

n. 

.39440  2 

o 

ni 

1.6528 

4 

1.26710 

e  =  1,0 

e  =  1/2 

,59510  3 

.4756  0  3 

1.57080 

1,0472  0 

1.2732  0 

1.1547  0 

,1488  0  2 

.049250 

1,3333 

.9455 

1.17840 

1.1075  0 

The  constants  for  e  =1.0  (the  circular  cone)  are  taken  from  Section 
6. IF.  These  results  will  be  studied  in  more  detail  to  determine  the  range  of 
Q  's  over  which  they  are  reasonably  valid.  The  pitching  moments  and  some  of 
the  damping  derivatives  for  the  elliptic  cone  will  also  be  evaluated. 

5.0  Analytic  Equations  for  the  Stability  Derivatives  Considering  Bluntness 

5.1  Introduction 

All  the  body  shapes  that  have  been  analyzed  up  to  this  point  have 
had  surfaces  whose  shape  can  be  expressed  by  a  single  equation.  It  has  also 
been  stated  elsewhere  in  this  report  that  the  methods  are  still  applicable  if 
the  top  and  bottom  of  the  body  can  each  be  expressed  by  different  equations, 
since  the  top  and  bottom  of  the  body  are  analyzed  separately.  These  general 
relationships  also  hold  for  bodies  whose  surfaces  are  composites  of  two  or 
more  shapes,  such  as  a  spherical  cap  and  a  cone.  The  only  details  to  be  studied 
are  how  the  results  for  the  different  shapes  are  to  be  combined  to  obtain  the 
aerodynamic  characteristics  of  the  complete  body. 

Surfaces  of  special  i uteres t  at  hypersonic  speeds  are  those  with 
blunted  noses  and  leading  edges.  Some  preliminary  investigations  have  been 
made  of  the  application  of  the  general  methods  shown  in  the  "First  Quarterly 
Technical  Report"  to  a  blunted  cone,  and  these  results  are  presented  in  this 
section. 


JO 
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5,2  Analytic  Equations  -  The  .Spherical  Cap 

In  order  to  analyze  the  aerodynamic  characteristics  and  stability 
derivatives  of  a  Spherically  blunted  cone ,  it  is  necessary  first  to  determine 
the  properties  of  the  spherical  cap.  Analytic  equations  for  the  force  and 
moment  coefficients,  and  the  stability  derivatives  of  a  spherical  cap  can  be 
determined  in  the  same  manner  as  that  used  for  the  cone  or  any  Other  surface. 
The  basic  methods  are  presented  in  Sections  3.0F,  4. OF,  and  5. OF,  The  basic 
equations  will  not  be  repeated,  only  the  values  of  the  constants  in  these 
equations  are  tabulated  below  for  the  bottom  half  of  a  cap  at  positive  angles 
of  attack. 


z 


Equation  of  the  Spherical  Cap  Segment 


-  {K  -  R)S-  y*  -  2*  *  R  =  3 


longitudinal  Force  and  lioment  Coefficient  Constants  for  C^ ,  C^.,  anti 


Longitudinal  Damping  Constants  for  C_ 


-*[(/-  -$-)  +z[t-  ]  +  !(/-  ^X%)(0-W)- 


^)  +  +  £  (/-  ty(2fc){Z6  *  tan  d)  to*  9  - 


“*(/-  9  +  |r(/-  %)(4)<9 


(K^)z~  if)  +  (^)  ]  jp*“ 


[*('-# ■)  *  (■*->] 

-  ^  tar? 9  +  tar? 9  -  £  0  ta?Q  X 


Lateral-Directional  Damping  Constants  for  CKr<  C|p_.  Cn  ,  and  C^r 


(*nr)f  =-£(*- jf-)  9 


[«'-  y(^)] 


■+■  4?  4a.a  $  —  ^  -t  ^an  ^  — 


t*n* Q 
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All  the  static  moment  constants  are  zero  since  the  center  of  moments 
is  taken  at  the  center  of  radius  of  the  spherical  cap..  The  simple  C.G.  trans¬ 
fer  relatic  ips  of  Section  3.2F  and  4.2F  can  he  used  to  obtain  moments  and 
moment  derivatives  about  any  other  C.G.  location.  The  force  and  moment  coef¬ 
ficient  constants  are  based  on  the  base  area  (  Ag)  and  base  radius  (RjJ  of  the 
bottom  half  of  the  spherical  cap  segment. 

It  is  interesting  to  note  that  for  a  C.G.  located  at  the  center  of 
radius  of  the  spherical  cap  (x0  =  R,  zq  =  0)  all  the  damping  constants  are 
zero  as  they  should  be.  Since  simple  C.G.  transfer  relationships  do  not  exist 
in  general  for  the  damping  derivatives,  the  constants  must  be  expressed  in 
terms  of  the  C.G.  position. 

For  the  case  of  a  half  sphere  ( &  -  0) ,  the  static  force  coefficient 
constants  assume  an  especially  simple  form 

(%a)i  *  (*V»)s  =  (*vl?  *  =  / 

=  (**)/  =  -jr 

(*Y p)t  ~  = 

Fi'om  these  constants,  it  is  possible  to  write  an  equation  for  C^, 

Cv,  and  Cv  of  the  bottom  of  the  half  spherical  cap  at  positive  angles  of 
attack. 


(5.1) 
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CT/3  =  " VST  s"> {*5 * * *)  (5 . 3) 

From  these  equations  it  is  easy  to  show  that  at  c*  =  o  ,  Cy  =  -C.,  ,  a 

jS  ^ 

not  unexpected  result. 

5.3  Longitudinal  Force  Coefficients  of  a  Blunt  Cone 

For  the  present,  only  equations  for  the  longitudinal  force  coeffi¬ 
cients  of  a  blunt  cone  have  been  derived.  The  moments  and  damping  derivatives, 
both  longitudinal  and  lateral-direction,  will  be  treated  in  the  next  quarter. 

In  many  respects  the  procedure  is  basically  that  which  will  be  shown  below. 

The  primary  additional  complication  is  that  static  moments  and  damping  deriva¬ 
tives  of  various  parts  of  a  composite  configuration  are  a  function  of  arm 
lengths  or  distance  from  the  center  of  moments  (C.G.)  of  the  total  configuration. 

In  order  to  derive  the  longitudinal  force  coefficients  of  a  blunt 
cone,  certain  area  relationships  are  required.  These  relationships  are  presented 
and  explained  with  an  appropriate  diagram  below. 


it1 
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The  bluntness  factor  (  ffi  )  of  the  cone  will  be  defined  in  terms  of 
the  base  radius  of  the  cone  (  )  and  the  base  radius  of  the  spherical  cap 

segment  (JL_  ). 


4  =  (4?)' 


The  total  planform  area  of  the  blunt  cone  (  S  )  is  composed  of  the  truncated 
cone  planform  area  (S^  )  and  the  planform  area  of  the  spherical  cap  (  ). 


s*  +  Se 


The  cap  force  coefficients  are  based  on  the  cap  base  area  (  /.  ,:j  ).  The  truncated 
cone  force  coefficients  are  based  on  the  truncated  cone  planform  area  (  ). 

If  the  blunt  cone  force  coefficients  tire  based  on  the  total  area  (  S  ),  then 
the  important  area  ratios,  expressed  in  terms  of  the  bluntness  factor  and  the 
cone  angle,  are  found  to  be  the  following: 

4? - J - f  (5.6) 

/  r  jjr 6  [/  -  ^  (0  ±  s/nZe) J 


Sc  . 

~r- 


/  *  )  (an  Q  secV  [/-  ^  (&  *  ^  SM  l&)  ~\ 


Based  on  total  planform  area,  the  axial  Force  coefficient  (C^)  and 
normal  force  coefficient  (C^)  of  the  blunt  cone  can  be  written  as 

Cfi  ~  (m)s  f— §e-)+  (cA)c  (-% ) 


<-,V  ~ (<-  m)&  (“X')  +  W  V  ( ~q~) 


r 
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The  subscripts  B  and  C  refer  to  the  spherical  cap  and  truncated  cone  respect¬ 
ively.  In  terms  of  the  constants  derived  for  the  cone  and  the  cap  these 
equations  become 

Q  =  [(**)/&  Cos*<*  ±{Ka)zb  *  (^4)38  &,no(  <***]  (^)  + 

■»“ ),c  coslx  +  {Ka)zc  $">*<*  +  [Ka)3C  s,n<*  Cos<*\{^t) 

\ 

Cn-  [(**)/£  C0S  01  *  (*A/)z8  SwZ<x  +  (Ayv)*3  Si  nek  <-<?$<*](  ^-)  + 

+[(*#)/ C  +[*A/)zc  5//7^  +  (Av)3c 

One  may  define  the  composite  or  blunt  cone  constants  as  follows: 


(*a)/GC  ~  +  (*A )/C  (~lP)  (5-8^ 

{*aI*£C  ~  (Ka)£s(ts’)  +  )j!C  (%)  (5.9) 

{<dWc  *  (*4 )y<T  (if)  (5.10) 

(%~)  (5.11) 


(<ft)zsc  +  (*v)/c  (%) 

IS 


(5.12) 
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(^CvJsBC  ~  {/<A/)3b{~s)  +  (5.13) 

The  arid  equations  now  become, 

J& 

Ot  =  (XA  )/BC  Cos  *  +  {XA)zac  stn  o(  +  {ka  jsgc  ease*  (5,14) 


2. 

c A/ -(**)/ ec  Cos  *  +  (x#)*ec  S/n  *  +  (*#h*e  s">«  ces<X  (5.15) 

These  equations  for  the  blunt  cone  obviour.v  take  the  same  form  as  the 
equations  for  the  sharp  cone,  or  any  other  aerodynamic  shape.  The  only  differ¬ 
ence  is  in  the  values  of  the  constants.  Eqs.  (5.14)  and  (5.15)  can  also  be  ex¬ 
pressed  in  the  alternate  and  more  useful  form  discussed  ^ n  Section  3.3F. 


CA  *  [Q-o)bc  *(a/)ec  s,n  [(<9a)sc  +  °<]  (5.16) 

C#  =  (n0 )BC  *  {n, )BC  *">*  [(<9r>)ec  +«]  (5.17) 


where 


f^/W  -  [(*4)*^  ~(^5*)/acJ  HX*hac 

(a-e)ec~  x  [(*>)/4c  +  fo)**e]  ~  X  (a/)oc 


(3a)ec  -  £ 


-i 


M 


4  isec 


2  -* _ (*4  )jac _ 


-19 
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(^i)bc  ±^§'iw)zsc  ~  [Xt/bec]  +  (** )sec 
{n6)BC  ~  &  [(**W  *  Wzec]  ~  2"  (ni)BC 


(& h)ec. 


/  .  -i 
zto* 


(*A/)30C 


(K*)tac~  (*#)/ec 


,  -/  {*n)38C 

M  S/n  p5~=m,=:^==fs= 

V  {[K/jizBC  ~  (**Wc]  +  (*/v)3dC 


It  can  easily  be  shown  that  the  truncated  cone  constants  based  on 
planform  area  are  the  same  as  the  constants  for  the  sharp  cone  and  can  be  ob¬ 
tained  from  Section  6. IF.  The  constants  for  the  spherical  cap  are  those  shown 
in  Section  5.2.  Obviously  it  is  possible  with  the  procedure  outlined  to  deter¬ 
mine  equations  for  CA  and  with  numerical  constants  for  any  blunt  half  cone 
at  positive  angles  of  attack  since  the  complete  cone  surface  will  "see"  the 
flow.  The  procedure  to  be  followed  for  the  blunt  full  cone  at  all  angles  of 
attack  is  the  same  as  that  for  sharp  cones  as  discussed  in  Section  6. IF. 

It  is  possible  to  obtain  reasonably  simple  analytic  expressions  for 
the  constants  in  Eqs.  (5.16)  and  (5.17)  if  it  is  assumed  that  the  half  cone 
angle  (  9  )  is  small,  and  the  bluntness  factor  (  fj  )  is  also  small.  Re¬ 

taining  only  up  to  third  order  terms  in  9f g,  and  their  products  we  have 


(a*)ac  "rg) 

-  f)®']-? -/)o] 
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(n/)ec  5  he[n-fe-  ’Jf-fl’  (5.22) 

(&)ec  =  -  #V]  +  *]  (5.23) 

The  initial  terms,  involving  only  the  angle  0  ,  are  the  values  of  the 
constants  for  the  sharp  cone.  The  terms  containing  f  D  are  the  contributions  of 

D 

bluntness.  These  equations  for  the  constants  can  be  compared  to  the  exact  values 
computed  numerically  using  Eq.  (5.8)  through  Eq.  (5.13).  The  exact  values  are 
plotted  as  a  function  of  the  cone  half  angle  (  Q  )  in  Figs.  13  and  14.  The  ap¬ 
proximate  expressions  for  the  constants  will  give  good  results  for  f^  < .16  and 

9  <  12°. 


It  should  be  possible  to  develop  similar  expressions  for  the  moments 
and  damping  derivatives  of  the  blunt  cone.  This  will  be  discussed  in  some  de¬ 
tail  in  a  future  report. 

6.0  Hypersonic  Pressure  Relationships  -  Deviations  from  Newtonian 

6.1  General  Discussion 

All  of  the  work  that  has  been  done  on  stability  derivatives  up  to 
this  point  has  been  based  on  the  simple  Newtonian  pressure  coefficient  rela¬ 
tionship  at  a  point  2  s in~ 8  ,  where  S  is  the  angle  between  the  free 

stream  flow  and  the  plane  tangent  to  the  surface  at  the  point.  Under  the  as¬ 
sumptions  of  Newtonian  flow,  this  equation  applies  only  as  long  as  $  >  0  • 
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For  negative  S' s  ,  the  pressure  coefficient  is  assumed  to  he  zero.  In  the 
strictest  sense,  the  Newtonian  pressure  relationship  applies  only  as  M  -*■  «e 
and  the  gas  constant  if  is  assumed  to  he  one.  In  a  practical  sense,  the  New¬ 
tonian  relationship  appears  generally  to  give  reasonable  results  for  the  total 
forces  as  long  as  M  &  >  4.0,  and  the  body  is  highly  three-dimensional,  such  as 
a  cone.  For  the  negative  pressure  region  of  the  body  |M  J  >  4.0*  the  negative 
pressures  are  less  than  5  percent  of  the  Newtonian  positive  pressures,  and 
theref ore  can  be  safely  neglected. 


An  alternative  to  the  Newtonian  constant  2  has  been  suggested  which 
gives  better  results  for  the  pressures  on  blunt  noses.  The  relationship  is  of 
the  form 


CP=  K  sin2 8  (6.1) 

Where  the  constant  K  is  taken  as  the  stagnation  pressure  coefficient  behind  a 
normal  shock. 

Any  pressure  relationship  of  this  form,  that  can  be  assumed  to  apply 
for  the  entire  body  at  all  values  of  1!  S7  leads  to  no  special  problems.  All  of 
that  which  has  been  done  on  the  stability  derivatives  assuming  the  Newtonian 
constant  2,  is  equally  valid  for  any  other  value  of  K.  The  only  requirement  is 
that  the  equations  for  force  and  moment  coefficients  and  stability  derivatives 
must  be  multiplied  by  the  ratio  K/2. 

It  has  been  established  that  the  pressure  coefficients  of  similar 
bodies  can  be  related  to  one  another  at  hypersonic  speeds  through  the  similarity 
parameter  11  $  .  But  for  different  body  shapes,  no  such  simple  relationship 
exists.  In  a  strict  sense,  the  constant  TC  in  Eq.  (6.1)  that  is  applicable  at 
any  point  on  the  body  is  not  only  a  function  of  the  local  value  of  M  fi ,  but  is 
also  a  more  complicated  function  of  the  general  body  geometry.  The  pressures 
behind  a  blunt  nose,  for  example,  are  a  function  of  the  bluntness  as  well  as 
the  local  value  of  '!£'  .  A  flat  plate  delta  wing  at  angle  of  attack  with  con¬ 
stant  free  stream  H  &  ,  has  a  pressure  coefficient  which  in  general  varies 
along  the  span. 
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Obviously  the  problem  of  computing  pressures  at  hypersonic  speeds  on 
surfaces,  and  integrating  the  results  to  obtain  the  stability  derivatives.,  is 
an  extremely  complex  one.  A  first  approach  to  this  problem  is  to  consider  those 
configurations  which  obviously  have  uniform  pressures  everywhere,  and  try  to  es¬ 
tablish  the  proper  pressure  relations  as  a  function  of  only  the  hypersonic  simi¬ 
larity  parameter  ( M  S  ).  One  such  configuration  is  the  circular  cone  at  zero 
angle  of  attack,  another  is  the  two-dimensional  flat  plate  at  angle  of  attack. 

For  the  present,  only  those  conditions  with  nose  and  leading  edge  shocks  attached 
Will  be  considered.  The  Newtonian  pressure  coefficient  relationship,  Eq.  fe.l) , 
and  the  combined  hypersonic-supersonic  similarity  law  suggest  that  the  form  the 
pressure  coefficient  might  take  is 


(6.2) 


where 


S  =  angle  between  the  surface,  and  free  stream 

(3  =  4mz-1 

For  sufficiently  high  Mach  numbers  and  small  S' s,  Eq.  (6.2)  will  reduce  to  the 
form 


Obviously  the  assumptions  of  Newtonian  flow,  V  1.0  and  y 
^(/Ssin  <J)  2. 


(6.3) 

— *■  oo  ,  are  such  that 


Eq.  (6,2)  has  some  interesting  implications  in  the  light  of  the  force 
and  moment  coefficient  and  stability  derivative  equations  that  have  been  derived 
assuming  Newtonian  flow.  All  of  these  results  have  shown  that  the  force  and  mo¬ 
ment  coefficients  are  a  Function  of  a  parameter  of  the  fora  sin*;  (  and 

the  derivatives  are  a  function  of  sin  (  The  angle  6^  can  toe  considered 
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to  be  an  average  integrated  surface  angle,  or  the  "effective"  surface  angle 
which  varies  with  the  surface  shape  and  the  force  coefficient  or  derivative 
being  evaluated.  In  the  light  of  Eq.  (6.2),  it  may  be  possible  to  consider 
that  for  the  entire  surface  being  analyzed,  the  average  S  is  (  &x  +  ct) .  If 
this  is  possible,  then  the  surface  average  pressure  law  that  is  applicable  is 
obtained  by  simple  substitution  in  Eq.  (6.2) 


If  evaluation  of  Eq.  (6.4)  results  in  a  number  different  than  the  Newtonian  2, 
the  implication  is  that  the  computed  Newtonian  results  for  the  force  and  moment 
coefficients  and  stability  derivatives  at  (  <9X  *■  )  should  be  multiplied  by  the 

following  ratio: 


?  [/3  4-*)] 

Z 

Of  course  this  procedure  is  all  predicated  on  the  fact  that  Eq.  (6.2)  is  appli¬ 
cable  on  the  average  to  the  shape  under  consideration. 

3rt  some  hypersonic  speeds  and  angles  of  attack,  negative  pressures  on 
a  body  are  far  from  negligible.  If  an  applicable  equation,  similar  to  Eq.  (6.2) 
can  be  derived  for  expansion  angles  (  8  e  K  then  it  may  be  possible  to  account 
for  the  negative  pressures  that  are  neglected  in  the  simple  Newtonian  theory. 
For  negative  pressures  Eq.  (6.2)  becomes 

U'-5) 

It  has  been  established  that  within  limits,  the  function  g'  (  (3 sin  Sv  )  can 
certainly  be  evaluated  for  rrandEL -Meyer  expansion  on  a  two-dimensional  flat 
plate  at  hypersonic  velocities.  Whether  similar  relationships  are  possible  on 
other  shapes  has  not  as  yet  been  determined. 
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Eq.  (6,5)  also  has  interesting  implications  in  light  of  what  has  been 
done  in  Newtonian  flow.  When  the  angle  t  +  oi.  )  for  the  bottom  surface  or  the 
angle  (  -  o(.)  for  the  top  surfaces  becomes  negative,  the  procedure  has  been  to 

drop  the  negative  angle  term  completely  in  determining  force  and  moment  coeffici¬ 
ents  and  stability  derivatives.  Eq.  (6.5)  establishes  a  possible  method  for  con¬ 
sidering  such  negative  pressures  following  a  procedure  identical  to  that  Used  for 
positive  pressures.  Thus  the  possibility  exists  that  the  Newtonian  results  can 
be  corrected,  at  least  approximately,  to  account  for  negative  pressures  oh  the 
basis  of  the  "effective1'  negative  angle  C)  of  the  surface,  where 

Be  - [&% +  °0  for  the  bottom  when  oC  <  ~  0-y 

Se  =  [By  -  <*)  for  the  top  when  c<  > 

The  procedures  suggested  above  for  correcting  the  Newtonian  stability 
derivatives  is  by  no  means  a  rigorous  one.  It  will  of  course  be  checked  by  com¬ 
paring  results  with  actual  experimental  data.  For  this  procedure  to  be  fruitful, 
appropriate  relationships  of  the  form  of  Eq.  (6.2)  and  (6.5)  must  be  found.  This 
approach  is  only  a  beginning,  and  as  the  program  progresses  other  methods  will 
also  be  investigated. 

6.2  Pressure  Laws  at  Hypersonic  Velocities 

The  present  discussion  will  be  concerned  with  methods  available  for 
establishing  the  functional  relationships  indicated  by  Eq,  (6.2)  and  Eq.  (6.5) 
at  hypersonic  velocities. 

For  the  two  dimensional  case-  oblique  shock  andPiaidtl  -Meyer  expansion, 
the  functional  relationships  of  Eqs.  (6,2)  and  (6-5),  with  certain  restrictions, 
have  been  established  and  appear  in  the  literature  in  various  analytic  forms. 

Tiider  the  assumption  that  the  shock  angle  and  wedge  angle  (  §  )  are 
small,  say  S  <  10l  t  and  M  >4.5-  it  is  possible  to  establish  a  very  simple 
expression  for  the  pressure  coefficient  of  the  wedge,  which  takes  the  form 
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Assuming  y  to  be  constant,  the  only  variable  on  the  right  hand  side  of  the 
equation  is  the  hypersonic  similarity  parameter,  M  £  ,  Eq.  (6.6)  can  be  ex¬ 
panded  in  various  power  series  that  give  good  results  for  particular  ranges  of 
H  <5  .  These  expressions  are  shown  below. 
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All  of  these  equations  give  results  with  less  than  about  2  percent  error  in  the 
ranges  specified,  and  can  be  checked  by  comparison  with  the  plot  of  Eq.  (6.6)  on 
Fig.  15.  Tt  is  obvious  from  this  plot  that  two-dimensional  pressures  are  much 
higher  than  those  predicted  by  the  simple  Newtonian  equation  at  the  lower  values 
of  21  $. 


Similar  expressions  can  be  established  for  lYanitl -Meyer  expansion  at 
hypersonic  speeds.  For  expansion  from  free  stream  Mach  Number  (M)  through  a 
negative  angle  .  the  pressure  coefficient  in  two-dimensional  flow  can  be 
expressed  as 
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The  only  important  assumption  required  for  Eq.  (6.10)  to  he  valid  is  that  the 
free  stream  Mach  Number  be  reasonably  large,  say  M  >  4.5.  Various  simplified 
expressions  can  be  derived  from  Eq..  (6.10)  that  are  valid  over  limited  ranges 
of  the  expansion  Similarity  parameter  (M  (5e).  These  expressions  are  listed 
below. 
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It  is  interesting  to  note  that  Eq.  (6.11)  differs  only  slightly  from 
Eq.  (6.7).  The  constant  in  the  last  term  is  about  10  percent  larger  for  a  Y 
of  1.4.  Eq.  (6.10)  is  plotted  as  Fig.  16,  and  the  accuracies  of  Eqs.  (6.11), 
(6.12),  and  (6.13)  can  be  determined  by  direct  comparison. 


Comparison  of  Figs.  15  and  16  indicates  that  at  small  values  of  H  6  , 
!1  S  -  .10,  expansion  and  compression  effects  are  equally  as  important 

in  determining  forces  and  moments  on  a  body.  At  H  8  =  |tt|r0|  l.O,  compression 
effects  are  approximately  three  times  as  great.  At  13  111  8C  ]  ■■  .5.0,  the 
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expansion  effects  are  negligible. 


The  relationships  discussed  so  far  are  two-dimensional,  and  most 
shapes  of  interest  in  hypersonic  flow  are  in  general  quite  three-dimensional* 

A  three-dimensional  shape  that  has  been  analyzed  in  some  detail  at  zero  angle 
of  attack  is  the  cone.  Cone  pressure  coefficients  computed  by  the  faylor-Maccoll 
theory  have  been  tabulated  by  Kopal  (Ref.  l) .  This  data,  for  cone  semi-angles 
of  less  than  10°  is  plotted  as  a  function  of  the  similarity  parameter  (38  in 
Fig.  15. 


Several  explicit  analytic  solutions  for  the  pressures  on  a  Cone  have 
been  determined  and  exist  in  the  literature.  A  simple  solution  for  a  slender 
cone  in  supersonic  flow  is 


=  ZJn  ~  '  (G.i-i) 

This  has  only  vei'y  limited  applicability.  A  comparison  of  Eq.  (6.14)  with  the 
results  of  Fig.  15  show  that  it  is  applicable  only  for  (3&  <  0.1.  A  second 
order  approximation  for  the  pressure  on  a  cone  at  zero  angle  of  attack  was  de¬ 
termined  by  Broderick  (Ref.  2). 
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For  small  angle  cones  (  8  <  10*1  ) ,  and  H  larger  than  2,  the  above  equation  can 
be  simplified  and  made  a  function  of  only  the  similarity  parameter  (  (38  ). 
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This  equation  can  be  compared  to  the  Kopal  results  of  Fig.  15.  The  Comparison 
is  good  for  |3  5  <  .50. 

Cole's  hypersonic  slender  body  theory  makes  it  possible  to  determine 
second  order  approximations  to  the  surface  pressures  on  slender  bodies.  The 
Second  order  approximation  is  essentially  a  Correction  to  the  first  order  approxi¬ 
mation,  which  consists  of  simple  Newtonian  plus  centrifugal  forces.  The  second 
order  hypersonic  solution  for  the  pressure  on  the  cone  at  zero  angle  of  attack  is 

CP  _  £/+ 3  +  / 

8*  Z(MSf  (6,17) 


Eq.  (6.17)  checks  quite  Well  With  the  Kopal  results  shown  in  Fig.  15  for 

ftb  >  i*5’ 


It  becomes  obvious  from  the  two-dimensional  and  cone  results  presented, 
that  geometry  has  an  especially  strong  influence  on  the  magnitude  of  the  pressures 
experienced  by  a  body  in  hypersonic  flow  at  low  values  of  the  similarity  parameter. 
At  ft  8  or  M  8  °f  0.1,  the  pressure  on  a  slender  cone  is  only  1/4  of  the  two- 
dimensional  value.  At  higher  values  of  M  S  ,  say  M  6  -  5.0,  the  cone  pres¬ 
sures  is  only  13  percent  less  than  the  two-dimensional  result.  It  is  also  true 
that  under  these  conditions  both  cone  and  two-dimensional  pressures  do  not  differ 
greatly  from  simple  Newtonian. 

The  cone  pressure  results  shown  here  will  be  used  to  determine  Mach 
Number  corrections  to  force  and  moment  coefficients  and  stability  derivatives 
for  the  cone  that  have  been  determined  assuming  simple  Newtonian  theory.  By 
comparing  these  results  to  experimental  data,  the  degree  of  validity  of  the  ap¬ 
proach  suggested  in  Section  6.1  will  be  established. 

An  investigation  of  other  body  shapes  is  obviously  required  in  order 
to  bridge  the  gap  between  the  pressures  on  a  two-dimensional  wedge  and  a  three- 
dimensional  cone.  But  second  order  solutions  on  other  than  bodies  of  revolution 
at  zero  angle  of  attack  are  difficult  to  obtain,  A  second  order  solution  has 
been  obtained  on  a  delta  planform  with  diamond  shape  cross-section  (Ref.  3), 
and  these  results  will  be  investigated  further. 


ad 
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The  discussion  up  to  this  point  has  been  concerned  with  corrections 
to  the  Newtonian  results  for  small  angles  of  attack  With  attached  shocks.  A 
similar  problem  exists  at  larger  angles  of  attack  where  the  shock  may  actually 
be  detached.  The  investigation  of  this  problem  will  also  be  undertaken, 

7.0  Projected  Work  Third  Quarter 

During  the  next  quarter,  many  of  the  results  shown  in  this  report  for 
the  force  and  moment  coefficients  and  stability  derivative  equations  in  Newtonian 
flow  will  be  extended.  Specifically  the  investigation  into  the  stability  deriva¬ 
tives  of  the  sharp  elliptic  cone  will  be  continued.  The  Work  on  the  stability 
derivatives  of  the  blunt  cone  will  also  be  extended.  As  explained  in  the  "First 
Quarterly  Technical  Report",  the  purpose  in  these  investigations  is  to  show  the 
application  of  the  general  methods  derived  rather  than  a  detailed  investigation 
of  specific  families  of  shapes. 

A  beginning  has  been  made  in  this  report  on  the  application  of  other 
hypersonic  theories  to  the  prediction  of  stability  derivatives.  This  investi¬ 
gation  will  continue.  In  Section  6.1  a  tentative  procedure  has  been  outlined  for 
extending  the  Newtonian  results  to  lower  hypersonic  Mach  numbers.  This  procedure 
will  be  applied  to  the  cone,  and  possibly  other  shapes,  and  the  results  compared 
to  existing  experimental  data. 

During  the  next  quarter,  a  concentrated  effort  will  be  made  to  develop 
methods  for  computing  the  stability  derivatives  under  trim  conditions,  a  primary 
goal  in  this  whole  hypersonic  stability  investigation.  The  basic  methods  and 
procedures  presented  in  this  and  the  previous  quarterly  report  lay  the  ground¬ 
work  for  this  investigation.  The  first  shape  to  be  investigated  under  trim  con¬ 
ditions  will  be  the  delta -planform  with  diamond-shape  cross-section,  since  this 
will  be  easiest  to  treat.  The  controls  to  be  used  for  trim  will  probably  be  a 
trailing  edge  deflected  flap.  To  the  degree  possible,  stability  derivatives 
with  trim  will  also  be  investigated  on  the  circular  and  elliptic  cone. 
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The  discussion  up  to  this  point  has  been  concerned  with  corrections 
to  the  Newtonian  results  for  small  angles  of  attack  with  attached  shocks .  A 
similar  problem  exists  at  larger  angles  of  attack  where  the  shock  may  actually 
be  detached.  The  investigation  of  this  problem  will  also  be  undertaken. 

7.0  Projected  Work  -  Third  Quarter 

During  the  next  quarter ,  many  of  the  results  shown  in  this  report  for 
the  force  and  moment  coefficients  and  stability  derivative  equations  in  Newtonian 
flow  will  be  extended.  Specifically  the  investigation  into  the  stability  deriva¬ 
tives  of  the  sharp  elliptic  cone  will  be  continued.  The  work  on  the  stability 
derivatives  of  the  blunt  cone  will  also  be  extended.  As  explained  in  the  "First 
Quarterly  Technical  Report",  the  purpose  in  these  investigations  is  to  show  the 
application  of  the  general  methods  derived  rather  than  a  detailed  investigation 
of  specific  families  of  shapes. 

A  beginning  has  been  made  in  this  report  on  the  application  of  other 
hypersonic  theories  to  the  prediction  of  stability  derivatives.  This  investi¬ 
gation  will  continue.  In  Section  6.1  a  tentative  procedure  has  been  outlined  for 
extending  the  Newtonian  results  to  lower  hypersonic  Mach  numbers.  This  procedure 
will  be  pnplied  to  the  cone,  and  possibly  other  shapes,  and  the  results  compared 
to  exis.JLng  experimental  data. 

During  the  next  quarter,  a  concentrated  effort  will  be  made  to  develop 
methods  for  computing  the  stability  derivatives  under  trim  conditions,  a  primary 
goal  in  this  whole  hypersonic  stability  investigation.  The  basic  methods  and 
procedures  presented  in  this  and  the  previous  quarterly’  report  lay’  the  ground¬ 
work  for  this  investigation.  Ihe  first  shape  to  be  investigated  under  trim  con¬ 
ditions  will  be  the  delta -plnnform  with  diamond-shape  cross-section,  since  this 
will  be  easiest  to  treat.  The  controls  to  be  used  for  trim  will  probably  be  a 
trailing  edge  deflected  flap.  To  the  degree  possible,  stability  derivatives 
with  trim  will  also  be  investigated  on  the  circular  and  elliptic  cone. 
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FIGURE  L  NEWTONIAN  CONSTANTS  FOR  A  HALF  CIRCULAR  CONE 
DAMPING  IN  PITCH  -  £6,  AT  <V»a)*0  >  Uo'«) « 0 

CONDITIONS'  CU  -9  *£  *  0,  0 


FIGURE  1  NEWTONIAN  CONSTANTS  FOR  A  HALF  CIRCULAR  CONE 
DAMPING  IN  FITCH  -  C.G,  AT  <Ve|*  */l  ,  tv*0)  VARIES 

CONDITIONS  0  £  *  0  ,  0 


HALF  CONE  ANGLE  -  8  -  DEGREES 

FIGURE  X  NEWTONIAN  CONSTANTS  FOR  A  HALF  CIRCULAR  CONE 
DAMPING  IN  PITCH  '  C..6  AT  (V*>*  I  ,  tV*o>  VARIES 

CONDITIONS  a*  -  0  ,  /9  *  0  ,  ^  *  0 

GJ 


IN  RADIANS 


FIGURE  A.  NEWTONIAN  CONSTANTS  FOR  A  HALF  CIRCULAR  CONE 
DAMPING  IN  PITCH  OF  A  SLENDER  CONE 

CONDITIONS:  ai-0,0*  0,^*0 

i>5 
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FIGURE  5.  NEWTONIAN  CONSTANTS  FOR  A  HALF  CIRCULAR  CONE 
DAMPING  IN  YAW-C  6,  AT  (x0/c)=0,  *0 

conditions:  a  >  -  e  ,  /a  =  o  ,  <£  =  o 

et> 


FIGURE  6.  NEWTONIAN  CONSTANTS  FOR  A  HALF  CIRCULAR  CONE 
DAMPING  IN  TAW  -  C.G.  AT  (^/c)  =  t/5  ,  %/t0)  VARIES 

CONDITIONS:  d>-e,/9*0,<£=0 
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HALF  GONE  ANGLE  ~8  ~  DEGREES 

FIGURE  7  NEWTONIAN  CONSTANTS  FOR  A  HALF  CIRCULAR  CONE 
DAMPING  IN  YAW  -  CG  AT  U<>/c)  =  I  ,  VARIES 

CONDITIONS  a>-«,  |9*0  ,  <f>*o 
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LOW-  MOMENT  CENTER  ~  U0/c) 

FIGURE  t.  NEWTONIAN  CONSTANTS  FOR  A  HALF  CIRCULAR  CONE 
DAMPING  IN  YAW  OF  A  SLENDER  COW 

CONDITIONS:  <2  £  -•  ,  JS  *  0  ,  *  0 
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HALF  COME  ANGLE  ~  6  ~  DEGREES 

FIGURE  3.  NEWTONIAN  OONSTANTS  FOR  A  HALF  CIRCULAR  CONE 
DAMPING  IN  ROLL  -  C,6,  AT  *  0 ,  bo*)  VARIES 

CONDITIONS:  ai*l,  /6»  0,  0 


TO 


C«p  =  mnp  *'n  t®np  * 

Q|r  =  m<f.  sin  {©|j  +a) 

Assumed  Moment  Center 

Long.  ~  (x0/c)  =  0 

Lot.  ~  (»0/c)  =  0 

■ 
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HALF  CONE  ANCLE  ~  6  ~  DEGREES 

FIGURE  10.  NEWTONIAN  CONSTANTS  FOR  A  HALF  CIRCULAR  CONE 
DAMPING  IN  YAW  DUE  TO  ROLLING  AND  ROLL  DUE  TO  YAWING 
06.  AT  Up/®*  «  0  ,  ly0/e>  *  0  ,  Up/fe)  VARIES 
CONDITIONS;  d£-6  ,(9*0  ,  $*0 

n 


(3/°*)  ,  (9/°g) 


Cnp=  mnp  Sin(0np  +  a) 
Cjf  =  m^f  sin(8|r  +  a ) 


Assumed  Moment  Center 
Long.  ~  (xp/c)  =  2/5 


HALF  CONE  ANGLE  ~  9  ~  DEGREES 


FIGURE  IL  NEWTONIAN  CONSTANTS  FOR  A  HALF  CIRCULAR  CONE 
DAMPING  IN  TAW  DUE  TO  ROLLING  AND  ROLL  OUE  TO  YAWING 
C«G,  AT  (VC)*  t/3  »  *  0  »  <VC>  VARIES 

CONDITIONS;  a>-9,/3*0,<£*0 


TiVc) 


HALF  CONE  AN6LE  ~  8  ~  DEGREES 
FIGURE  \Z ,  NEWTONIAN  CONSTANTS  FOR  A  HALF  CIRCULAR  CONE 
DAMPING  IN  YAW  DUE  TO  ROLLING  AND  ROLL  DUE  TO  YAWING 
C  G.  AT  (stp/c) *  I  ,  (»0/c)  *  0  ,  VARIES 

CONDITIONS'  <X>~8,  j3  =  0,<£*0 
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HALF  CONE  ANGLE  ~6~  DEGREES 

FIGURE  13.  NEWTONIAN  CONSTANTS  FOR  A  BLUNT  CIRCULAR  HALF  CONE 

AXIAL  FORCE  COEFFICIENT 
CONDITIONS;  <I£-G,  /S*0  ,£'0 
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HALF  CONE  ANGLE  ~  9  ~  DEGREES 

FIGURE  14.  NEWTONIAN  CONSTANTS  FOR  A  BLUNT  CIRCULAR  HALF  CONE 

NORMAL  FORCE  COEFFICIENT 

CONDITIONS;  di-G  ,/3*0  ,  <£*0 
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